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2009-2010  Annual  Progress  Report 

This  report  presents  the  specific  aims  and  accomplishments  of  our  prostate  cancer  research 
project  during  the  year  of  funding  sponsored  by  the  US  Department  of  the  Army.  It  covers  our 
activities  from  May  1 , 2009  to  August  31 , 201 0. 

1.  Introduction 

It  is  known  that  effective  treatments  for  prostate  cancer  are  highly  associated  with  accurate  and 
early  detection  of  prostate  cancer.  The  current  tools  for  detection  of  prostate  cancer  have  serious 
shortcomings.  So,  the  final  diagnosis  for  prostate  cancer  relies  on  needle  biopsy  through 
transrectal-ultrasound  guidance.  Since  ultrasound  is  not  able  to  show  any  cancerous  lesions,  the 
current  diagnostic  needle  biopsy  for  prostate  cancer  is  rather  “Blind”:  the  biopsy  samples  were 
collected  almost  blindly  without  knowing  whether  or  not  the  biopsy  lesions  are  highly  specious  for 
cancer.  The  current  diagnostic  method  in  needle  biopsy  has  three  major  deficiencies:  (1)  it  may 
miss  cancers  that  are  present  and  have  large  high-grade  lesion,  (2)  it  can  underestimate  the 
severity  of  cancer  that  present,  and  (3)  it  may  find  only  tiny  low-grade  cancers  that  likely  do  not 
cause  problems.  Because  of  these  deficiencies,  the  current  prostate  cancer  patients  are  either 
under  or  over  treated:  under  treatments  cause  patients’  lives;  over  treatments  cause  patients’ 
unnecessary  impotence  or  incontinence. 

The  hypothesis  for  this  study  is  that  a  multi-modal  optical  spectroscopic  method  and  an 
integrated  needle  probe  can  be  developed  for  guiding  needle  biopsy  for  prostate  cancer  diagnosis. 
Multi-modal  optical  measurements  to  be  utilized  for  the  study  are  (1)  light  reflectance  spectroscopy 
(LRS),  (2)  auto-fluorescence  spectroscopy  (AFS),  and  (2)  auto-fluorescence  life-time 
measurements  (AFLM). 

The  project  has  four  specific  aims: 

Aim  1_:  to  develop  a  multi-modal,  optical  spectroscopic  instrument,  which  allows  the 
measurements  of  (1)  LRS,  (2)  AFS,  and  (3)  AFLM.  The  proposed  system  will  be  portable,  can 
be  used  for  in  vivo  measurements,  and  collect  and  present  the  data  in  real  time. 

Aim  2:  to  integrate  the  optical  fibers,  which  collect  light  scattering  and  auto-fluorescence  from 
the  prostate  tissue,  into  a  transrectal-ultrasound,  needle-biopsy  probe.  In  the  development 
phase,  the  optical  signatures  of  prostate  cancer  can  be  collected  with  the  biopsy  tissues  and 
identified  along  every  tract  of  needle  biopsies. 

A[m  3:  to  collect  optical  signals  of  control  and  cancer  tissues  ex  vivo  along  with  the  regular 
human  needle  biopsy,  followed  by  classification  algorithm  development  to  discriminate  cancer 
tissues. 

A[m  4:  to  perform  in  vivo  measurement  from  human  subjects  to  obtain  the  accuracy  and 
sensitivity  of  the  integrated  probe  in  order  to  provide  real-time,  on-site,  improved  guidance  for 
prostate  cancer  tissue  biopsy. 

This  report  will  summarize  the  work  we  performed  from  May  1,  2009  to  August  2010  and  the 
corresponding  achievements  obtained  during  this  period  of  time,  as  given  below. 
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2.  Body  of  the  Report 

The  PI  and  her  research  team  have  made  significant  efforts  from  2009-2010  to  accomplish  the 
proposed  aims,  resulting  in  several  published  journal  and  conference  papers.  We  will  give  a 
brief  report  in  the  following  sub-sections. 

2.1  Background 

Prostate  cancer  is  the  most  commonly  found  male  cancer  in  the  United  States  and  is  the 
second  leading  cause  of  death  from  cancer  in  men  [1],  It  was  estimated  that  in  2009,  the 
number  of  new  prostate  cancer  cases  reported  would  be  192,280,  resulting  in  27,360  deaths  [2]. 
At  present,  digital  rectal  examination  (DRE),  prostate  specific  antigen  (PSA)  blood  test  and 
transrectal  ultrasound  (TRUS)  guided  biopsy  are  the  clinically  available  techniques  for  prostate 
cancer  screening  and  diagnosis.  Among  these,  current  gold  standard  for  prostate  cancer 
detection  is  TRUS  guided  needle  biopsy,  which  involves  resection  of  a  core  of  prostate  tissue 
with  the  guidance  of  an  ultrasound  probe.  However,  though  ultrasound  is  able  to  provide 
anatomical  images  with  high  spatial  resolution,  it  lacks  the  sensitivity  in  differentiating  tumor 
from  normal  tissue,  especially  at  early  stages.  This  deficiency  makes  TRUS  guided  biopsy  a 
rather  “blind”  procedure  involving  quasi-random  sampling  of  the  prostate  tissue  [3].  The 
implications  of  this  drawback  are  false  negatives,  and  oversampling  of  tissue  due  to  non-specific 
targeting.  A  standard  biopsy  procedure  involves  resection  of  10-12  cores,  and  it  was  believed 
that  by  increasing  the  number  of  core  biopsies,  the  detection  rate  could  be  increased  [4], 
However,  recent  studies  indicate  that  there  is  no  significant  improvement  in  detection  rate  by 
increasing  number  of  biopsy  cores  (which  is  referred  to  as  saturation  biopsy  [5,6].  Moreover, 
increasing  the  number  of  cores  increases  complications  associated  with  the  biopsy  procedure 
[4],  It  is  therefore  imperative  to  develop  techniques  that  can  identify  the  tumor  in-vivo  with  high 
specificity  and  sensitivity  and  improve  the  accuracy  of  biopsy.  It  is  even  more  crucial  to  develop 
a  cancer-targeted  diagnostic  technique,  which  can  identify  the  tumor  on-site  with  high  sensitivity 
and  specificity. 

Optical  spectroscopic  techniques  can  be  a  candidate  for  such  needed  tool  to  identify  cancer 
from  control.  They  may  be  implemented  through  needle  probe  geometry  and  will  be  extremely 
helpful  in  improving  the  diagnostic  procedure.  The  first  step  in  this  study  is  to  utilize  an  optical 
fiber  which  can  be  used  to  obtain  a  priori  information  by  mapping  the  prostate  tissue  first, 
followed  by  subsequent  biopsy  of  selected  lesions.  In  this  way,  we  can  increase  the  accuracy  of 
the  biopsy  procedure,  substantially  reduce  the  number  of  biopsy  cores,  and  thus  reduce  the 
amount  of  tissue  resected  as  well  as  minimize  the  associated  complications.  Specifically,  two 
techniques  have  been  implemented  in  such  geometry  in  the  past  year:  (1)  light  reflectance 
spectroscopy  (LRS)  and  (2)  auto-fluorescence  lifetime  measurement  (AFLM). 

LRS  has  been  previously  utilized  for  a  variety  of  clinical  applications  by  characterizing  different 
tissue  types  using  either  empirical  or  analytical  methods  [7,8,9,10,11,12],  Recently,  updated 
methods  have  been  developed  to  quantify  the  absolute  values  of  optical  parameters  from  LRS 
with  small  probe  geometries  [9,13].  These  analytical  approaches  make  it  more  robust  to  be 
applied  as  a  generalized  means.  On  the  other  hand,  various  research  groups  have  also 
implemented  laser-  induced,  steady-state  auto-fluorescence  spectroscopy  and  time-resolved 
lifetime  measurement  for  cancer  detection  [14],  Steady-state  fluorescence  measurements  are 
easy  to  perform,  however,  they  are  intensity-dependent  and  thus  sensitive  to  optical  properties 
of  tissue,  ambient  light,  and  excitation  light  intensity.  We  opted  for  temporally-resolved 
fluorescence  lifetime  measurements  as  they  are  independent  of  absolute  intensity  of  excitation 
light  and  fluorophore  concentration. 


4 


2.2  Report  for  Aim  1 


Aim  1:  develop  a  multi-modal,  optical  spectroscopic  instrument,  which  allows  the 
measurements  of  (1)  LRS,  (2)  AFS,  and  (3)  AFLM.  The  proposed  system  will  be  portable,  can 
be  used  for  in  vivo  measurements,  and  collect  and  present  the  data  in  real  time. 

For  this  aim,  we  have  implemented  and  evaluated  both  LRS  and  AFLM  systems  independently, 
to  differentiate  between  cancer  and  non-cancerous  tissue  using  laboratory  tissue  phantoms. 
LRS  is  used  to  calculate  the  absolute  concentrations  of  oxy-hemoglobin,  deoxy-hemoglobin  and 
scattering  properties  of  the  tissue.  AFLM  uses  a  custom  built  device  to  measure  auto¬ 
fluorescence  lifetime  of  the  tissue.  We  have  tested  the  efficacy  of  using  these  multiple 
parameters  as  classifiers  to  identify  prostate  cancer.  The  results  are  presented  in  this  report  and 
provide  us  with  the  footprints  for  future  ex-vivo  and  in-vivo  human  prostate  cancer 
measurements. 

2.2.1  Laboratory  tissue  phantoms  and  preparation 

Phantom  measurements  were  carried  out  to  test  our  newly  implemented  LRS  and  AFLM  system 
for  this  part  of  study.  Figure  1  shows  the  diagram  of  the  combined  system  with  a  tissue 
phantom.  To  create  a  homogenous  phantom,  a  blood-intralipid  solution  was  made  by  diluting  a 
stock  intralipid  solution  (concentration=20%,  Baxter  Healthcare  Corporation,  Deerfield,  IL)  with 
phosphate  buffered  saline  (PBS)  and  de-fibrinated  horse  blood  (Hemostat  Laboratories,  Dixon, 
CA).  We  varied  the  concentrations  of  intralipid  and  blood,  to  obtain  various  combinations  of 
scattering  and  absorption  properties.  To  simulate  hemodynamic  changes  in  the  phantom,  4  mg 
of  yeast  was  added  into  a  3-liter  liquid  phantom  to  deoxygenate  the  blood-containing  phantom. 
When  the  phantom  was  fully  deoxygenated,  100%  oxygen  was  slowly  bubbled  into  the  solution 
to  re-oxygenate  it.  Overall,  we  could  simulate  changes  in  light  scattering  by  alternating  intralipid 
concentrations  and  mimic  changes  in  light  absorption  by  modifying  either  the  blood  volume  or 
oxygenation  state.  In  the  tumor-simulating  experiments,  we  increased  the  light  absorption  by  2-3 
times  higher  than  that  from  the  baseline  (healthy-tissue)  phantoms  to  detect  the  differences  in 
both  LRS  and  AFL. 


Fig.  1  A  block  diagram 
illustrating  experimental 
set-up  for  LRS  system 
(left)  and  AFLM  system 
(right).  Measurements 
were  made  sequentially, 
by  placing  the  fiber  tips  on 
the  tissue  phantom 
surface.  A  closer  view  of 
bi-furcated  fiber  tips  are 
shown  for  each  fiber  (red  = 
source;  blue  =  detector). 


Tissue 

phantom 
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2.2.2  Auto-florescence  lifetime  measurement 


AFLM  Instrumentation 

For  auto-florescence  lifetime  measurements,  a  custom-made,  single-channel,  Time  Correlated 
Single  Photon  Counting  (TCSPC)  system  (ISS  Inc.,  Champaign,  IL)  was  employed.  Fig.  1  (right) 
shows  an  overall  experimental  setup  for  the  time  domain  AFLM.  The  system  consisted  of  a  5- 
slot  filter  wheel  for  excitation  wavelength  selection  (which  was  driven  by  a  12-V-powered 
stepper  motor),  an  emission  filter,  a  continuously  variable  neutral  density  (ND)  filter  for 
excitation  light  intensity  control,  and  a  highly  sensitive  cooled  PMT  (Becker  &  Hickl  GmbH)  with 
wavelength  sensitivity  between  125  -  850  nm.  The  PMT  gain  was  controlled  via  PC-based  card 
(DCC-100).  A  picoseconds  (ps)  broadband  (400  -  1800  nm)  pulsed  laser  (SC  -  450,  Fianium 
Inc.,  Eugene,  Oregon)  was  used  as  an  excitation  source  at  the  repetition  rate  of  20  MHz. 
Precise  synchronization  between  the  incoming  laser  pulse  and  photon  event  was  achieved 
through  a  PC-based  single  photon  counting  card  (SPC-130).  As  shown  in  Fig.  1,  the  laser  was 
coupled  to  the  source  fiber  of  the  bi-furcated  optical  probe  (core  diameter  100  pm);  resulting 
fluorescence  emission  was  collected  through  the  detector  fiber  (400  pm  core  diameter). 

AFLM  Data  Analysis 

In  this  part  of  study,  we  were  mainly  targeting  auto-fluorescence  from  Flavins,  Lipo-pigments 
and  Porphyrins  as  a  contrast  parameter  to  differentiate  cancer  from  healthy  tissue.  In  order  to 
achieve  maximum  possible  auto-fluorescence  from  all  three  endogenous  compounds,  an 
excitation  wavelength  of  447  nm  was  chosen  [15].  While  keeping  the  excitation  wavelength 
constant,  the  emission  wavelengths  were  changed  between  532  nm,  562  nm,  632  nm  and  684 
nm.  Data  collection  was  done  through  Vinci  software  provided  by  ISS.  The  optical  probe  was 
placed  at  ~1  mm  away  from  the  surface  of  the  tissue  phantom;  the  fluorescence  data  were 
collected  at  5  random  positions  from  the  phantom  for  each  emission  wavelength.  The  detected 
fluorescence  decay  was  imported  into  Matlab  (The  Mathworks  Inc.,  Natick,  MA)  and  normalized 
with  respect  to  peak  intensity. 

In  order  to  quantitatively  differentiate  auto-fluorescence  decay  of  cancer  and  control,  each  curve 
can  be  fitted  (least  square  non-linear  fitting  algorithm)  using  a  two-component  exponential 
model,  as  given  by: 


t 

Intensity(I)  =  Y.oCje  T‘  , 
i 


(1) 


where,  n{ i=1,  2)  indicates  lifetime  of  each  component  and  or,  (i=1,  2)  is  intensity  contribution  of 
each  component  to  the  overall  fluorescence  decay.  The  confidence  interval  for  each  of  these 
four  fitted  parameters  was  computed  using  the  “confint”  function  in  Matlab.  In  order  to  make  the 
comparison  between  normal  and  cancerous  tissue,  the  mean  lifetime,  <x>,  was  calculated  using 
equation  2: 


Z  aiXTi 


(2) 
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2.2.3  Light  Reflectance  spectroscopy 


LRS  Instrymentation 

The  LRS  system  consisted  of  a  tungsten-halogen  light  source  (HL2000HP,  Ocean  Optics, 
Dunedin,  FL,  USA),  a  CCD  array  spectrometer  (USB  2000+,  Ocean  Optics,  Dunedin,  FL  USA), 
and  a  laptop  computer.  The  spectrometer  provided  a  spectrum  ranging  from  460  nm-1150  nm. 
A  custom-made,  bi-furcated,  fiber  optic  probe  (Fiberguide  Industries,  Stirling,  NJ,  USA)  was 
designed  and  implemented  with  source  and  detector  diameters  of  100  pm,  center-to-center 
source  and  detector  separation  of  -100  pm,  and  an  outer  probe  diameter  of  800  pm.  The  probe 
was  fixed  on  a  stereotactic  frame  in  order  to  control  the  probe  tip  position  and  minimize  the 
pressure  on  the  tissue  surface.  Multiple  data  points  were  obtained  by  placing  the  probe  at 
different  spatial  locations  on  the  healthy-tissue  phantoms  and  then  the  tumor-simulating  tissue 
phantoms.  On  an  average,  5  data  points  were  obtained  for  each  tissue  phantom  type. 


LRS  Data  Acguisition  and  Analysis 


Each  acquired  spectrum  was  divided  by  a 
reflectance  spectrum  obtained  from  a  diffuse 
reflectance  standard  (WS-1,  Ocean  Optics, 
Dunedin,  FL,  USA)  to  eliminate  the  spectral  effect 
from  the  light  source,  optical  fibers  and  CCD 
array  detector  from  the  reflectance  spectrum 
taken  from  the  tissue  sample.  Figure  2  shows  a 
schematic  diagram  for  the  setup.  For  system 
calibration,  we  utilized  a  gold-standard,  dual¬ 
channel,  tissue  oximeter  (ISS,  Inc.,  IL,  USA),  as 
also  shown  in  Fig.  2.  The  ISS  oximeter  could 
provide  absolute  values  of  optical  properties  for 
the  measured  sample  under  light  interrogation, 
thus  providing  calibration  parameters  to  LRS. 


Then,  the  spectral  range  from  500  nm  -  850  nm  was  selected  and  a  reflectance  model  [9]  as 
shown  in  equation  (3)  was  utilized  to  obtain  absolute  values  of  the  concentrations  of  oxy¬ 
hemoglobin  ([HbO]),  deoxy-hemoglobin  ([HbR])  and  total  hemoglobin  ([HbT]  =  [HbR]  +  [HbO]), 
along  with  the  reduced  scattering  coefficient  ( fjs ). 


RP(X)  = 


k\  +k2  -AaW  ’ 


(3) 


where  Rp  is  the  measured  reflectance,  fjs’  is  the  reduced  scattering  coefficient,  pa  is  the 
absorption  coefficient  depending  on  the  concentration  of  hemoglobin  derivatives  and  /q  and  k2 
are  calibration  constants  that  depend  on  the  probe  geometry  and  hardware  set-up.  A  detailed 
description  of  methods  to  obtain  the  above  parameters  has  been  previously  published  [7]. 
Briefly,  at  first,  /q  and  k2  were  calculated  for  the  specific  fiber  optic  probe  using  tissue  phantom 
calibration.  Then,  an  optimization  routine  [16]  (ant  colony  optimization)  was  used  to  obtain  the 
parameters  of  /js’  at  750  nm,  [HbO]  and  [HbR],  by  achieving  the  best  fit  for  the  measured 
reflectance.  The  data  processing  was  implemented  and  executed  in  Matlab.  Specifically,  we 
provide  details  on  data  analysis/process  procedures  as  follows. 


(1)  Quantification  of_  hemoglobin  concentrations : 


7 


It  is  well  known  that  /ja(A)  is  a  function  of  concentrations  of  deoxygenated  hemoglobin,  [Hb], 
oxygenated  hemoglobin,  [HbO],  and  water,  H20.  The  spectral  dependence  of  ga  on  [HbO],  [Hb], 
and  H20  for  blood-perfused  tissues  can  be  written  as 

Ma{X)  =  [HbO]  sHboU)  +  [Hb]  sHh(X)  +  [%H20]sHi0a),  (4) 

where  ‘X  is  the  wavelength  in  nm,  eHb0  (A),  eHb  W  and  eH20  (A)  are  the  wavelength  dependent 
extinction  coefficients  of  [HbO],  [Hb],  and  water,  respectively.  [%H20]  represents  the  percentage 
of  water  in  the  medium.  Other  chromophores,  such  as  fat  and  melanin,  may  be  added  to  this 
equation,  depending  on  the  wavelength  range  of  interest  and  tissue  type  under  investigation. 
The  spectral  dependence  of  gs  can  be  approximated  as  given  by  Mie  theory  [17]: 


g  =1.1-(0.58x10_3)A, 

(5) 

Ms  =Ms0-~8), 

(6) 

£ 

II 

a 

<5- 

(7) 

where  gs  is  the  effective  scattering  coefficient,  gs’  is  the  reduced  scattering  coefficient,  and  g  is 
the  anisotropy  factor.  In  eq.  (7),  parameters  of  a  and  b  are  constants  and  depend  on  scatterer 
sizes  and  types.  For  10%  intralipid  solution,  the  calculated  values  are  a  =  2.54x1 09  cm"1  and  b  = 
2.4  [17]. 

(2)  Calculation  and  calibration  of  k1  and  k9 

In  this  study,  we  determined  the  values  of  /q  and  k2  using  the  least-squares  regression 
approach.  By  rearranging  (3),  we  obtain 

=  kx+  k2Ha (T)  ■  (8) 

KP  (A) 

This  equation  shows  that  /q  and  k2  can  be  determined  by  obtaining  a  linear  regression  line 
that  best  fits  gs’(A)/Rp(A)  versus  ga  (A).  In  principle,  we  can  obtain  the  measured  spectra  of  RP(A) 
from  the  LRS  system  within  500-850  nm  and  also  achieve  quantification  of  ga  and  gs’  at  750  nm 
and  830  nm  from  ISS  oximeter  if  the  measurements  by  both  LRS  and  ISS  oximeter  are  taken 
simultaneously.  Besides  ga  values,  ISS  oximeter  is  also  able  to  provide  derived  values  of  [Hb], 
[HbO],  and  gs’  at  two  wavelengths  (OxiplexTS).  Then,  wavelength-dependent  absorption 
spectra,  / ua(A),  can  be  quantified  using  eq.  (4)  if  the  measured  [Hb]  and  [HbO]  as  well  as  the 
hemoglobin  extinction  coefficients  are  available.  Moreover,  given  two  gs’  values  at  two 
measured  wavelengths  (750  nm  and  830  nm),  it  is  reasonable  to  interpolate  and  extrapolate 
/ us’(A)  over  the  desired  wavelength  range  using  eqs.  (5)-(7)  based  on  Mie  theory. 

During  the  system  calibration  in  the  study,  we  followed  the  exact  principle  or  procedures 
given  above  to  acquire  RP(A),  ga(A),  and  gs’(A)  from  a  set  of  blood-containing  tissue  phantoms.  A 
corresponding  set  of  /q  and  k2  at  various  total  hemoglobin  concentrations  and  reduced 
scattering  coefficients  were  obtained,  averaged,  and  calculated  for  their  means  and  standard 
deviations.  If  the  relative  errors  for  both  /q  and  /qwere  less  than  10%  (i.e.,  standard  deviations 
divided  by  their  means)  ,  then  the  corresponding  set  of  k1  and  k2  were  finalized  as  the  system 
calibration  parameters  for  the  specific  LRS  system  (i.e.,  the  light  source,  fiber  probe,  and 
spectrometer)  chosen  for  this  study. 

As  an  example,  Fig.  3(a)  shows  a  set  of  measured  RP(A),  ga  (A),  and  gs’(A)  (blue  symbols) 
and  the  linear  regression  line  (red  thick  line)  for  our  LRS  probe  with  /q  being  the  y-intercept  and 
k2  being  the  slope.  While  the  inverse  calculations  for  an  unknown  sample  will  be  described  in 
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the  next  sub-section,  an  illustration  of  the  measured  and  fitted  reflectance  is  shown  in  Fig.  3(b). 
The  accuracy  of  parameters  /q  and  k2  critically  affects  the  ‘goodness  of  fit’  in  the  inverse 
calculations,  and  thus  influences  the  accuracy  of  measured  optical  properties  of  tissues  under 
interrogation.  It  is  noteworthy  to  point  out  that  k1  and  k2  depend  highly  on  the  white  sample 
measurement  for  eliminating  instrumentation  effects  on  the  measured  spectrum  of  the  sample 
under  examination.  Incorrect  measurement  of  the  white  sample  will  lead  to  serious  errors  for  the 
fitted  parameters,  i.e.,  [Hb],  [HbO],  a  and  bin  eq.  (7). 


Figure  3 


White  sample  measurement  can  vary  with  the  following  parameters:  distance  between  the 
fiber  tip  and  sample  surface,  integration  time,  light  source  intensity.  Measuring  different  types  of 
tissue  requires  a  different  integration  time  for  the  detector.  Fig  4 
shows  the  measurement  geometry  from  the  white  sample.  Since 
the  reflectance  curve  needs  to  be  divided  by  the  reference  curve 
from  the  white  sample,  care  should  be  taken  to  standardize  the 
procedure.  This  can  be  done  either  by  fixing  the  distance  and 
collecting  white  sample  at  various  integration  times  to  create  a  look¬ 
up  table,  or  by  scaling  the  white  sample  reflectance  curve  based  on 
one  single  measurement.  The  latter  method  is  based  on  our 
observation  that  the  relationship  between  reflectance  spectra  and 
integration  times  is  approximately  linear. 

(3)  Inverse  calculation 

The  inverse  calculation  of  eq.  (3)  deals  with  determination  of  tissue  chromophore  concentrations 
and  effective  scattering  coefficients  from  the  measurements  of  LRS  taken  on  the  tissue  surface. 
The  mathematical  problem  of  inversion  can  be  practically  solved  by  an  approach  of  function 
minimization  using  an  optimization  algorithm.  The  algorithm  searches  for  an  optimal  set  of 
values,  as  represented  by  x=(  [Hb],  [HbO],  a,  b),  from  the  given  multi-parameter  space  that  best 
match  the  measured  Rp() 1)  using  the  least-  squares  analysis.  It  can  be  mathematically 
expressed  as  follows: 

f(x)  =  £ [Rp iXfi~ed)  - Rp (Aiypredicted) \  ,  (9) 

1=1 

where  f(x)  represents  an  objective  function  and  needs  to  be  minimized  by  finding  an  optimal  set 
of  the  parameters,  xopU  M  is  the  number  of  wavelengths,  Rp  (x$measured)  and  Rp  (Xi)(predicted>  are  the 
reflectance  values  (at  wavelength  A,)  measured  by  the  optical  probe  and  predicted  by  eqs.  (3)  to 
(7),  respectively,  with  the  parameter  set,  x,  to  be  optimized.  Evolutionary  algorithms  are  widely 
utilized  in  the  field  of  optimization.  The  algorithms  work  on  a  population  of  values  rather  an  initial 
guess,  making  the  solution  independent  of  an  initial  guess.  Their  inherent  characteristics  make 
them  suitable  to  search  for  a  global  minimum,  rather  than  being  stuck  in  local  minima.  The  ant 
colony  optimization  (ACO)  algorithm,  introduced  by  M.  Dorigo  [18],  is  a  probabilistic  evolutionary 
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Figure  4 


9 


technique  for  solving  computational  problems.  This  basic  technique  was  modified  in  this  study  to 
suit  function  minimization  [16]. 

Equation  (3)  is  very  non-linear  and  spanned  from  500  nm  to  850  nm.  In  order  to  obtain 
optimal  fit  between  the  measured  and  calculated  reflectance  curve,  we  had  to  incorporate 
multiple  sequences,  within  each  of  which  the  objective  function  for  minimization,  eq.  (9),  was 
selected  somewhat  differently.  First,  the  entire  spectrum  (500  nm  -  850  nm)  was  utilized  for 
fitting  and  the  fitted  values  of  [Hb],  [HbO],  a  and  b  in  eq.  (7)  were  obtained.  Next,  the  [js(A) 
values  (or  a  and  b)  were  fixed  with  15%  bounds  used  in  sequence  one.  A  smaller  spectral 
region  of  520-590  nm  that  has  a  strong  hemoglobin  absorption  band  was  selected  for  re-fitting 
to  improve  ‘goodness  of  fitting’  for  [Hb]  and  [HbO].  Finally,  after  refining  the  fitted  parameters  of 
[Hb]  and  [HbO],  they  were  fixed  within  20%  bounds  used  in  sequence  two.  The  entire  spectrum 
(500-850  nm)  was  re-fitted,  leading  to  the  finally  optimized  values  of  [Hb],  [HbO],  a  and  b. 
Overall,  three  fitting  sequences  with  changing  parameter  bounds  and  spectral  regions  were 
utilized  in  order  to  obtain  an  optimal  fit  for  scattering  and  absorption  parameters.  An  example  to 
show  comparison  between  the  measured  and  fitted  reflectance  is  already  given  in  Fig.  3(b). 

2.2.4  Results: 

Algorithm  validation  for  absolute  guantjfication  by_  LBS 

In  order  to  validate  the  algorithm  for  absolute  quantification,  we  created  tissue  phantoms  and 
performed  the  measurements  with  both  our  single-channel  LRS  system  and  a  tissue  oximeter 
(ISS  Oximeter),  which  was  used  as  a  gold  standard  device  for  comparison.  Figure  5(a)  shows 
comparisons  between  the  values  obtained  from  ISS  tissue  oximeter  and  those  quantified  with 
our  method  when  a  change  in  oxygen  saturation(02  Sat)  was  achieved  by  adding  yeast  to  the 
liquid  phantom.  The  measurements  were  made  at  12  different  time  points  in  the  02Sat  range  of 
~20%-80%  during  the  oxygenation  and  deoxygenation  process.  It  is  clear  from  Fig.  5(a)  that 
oxygen  saturation(02  Sat)  values  measured  by  both  modalities  are  very  consistant  in  the  02  Sat 
range  of  20%-80%.  The  absolute  errors  were  also  calculated  at  each  data  point  (n=12)  and  are 
plotted  in  Fig.  5(b):  average  absolute  errors  for  02  Sat,  [Hb]  and  [HbO],  along  with  standard 
error  of  mean  (SEM)  were  3.0±0.7  %,  0.9±0.2  pM,  and  0.9±0.2  pM,  respectively.  Similar 
experiments  were  repeated  several  times,  and  the  results  were  very  similar  to  those  seen  in 
Figs.  5(a)  and  5(b). 

The  reliability  of  the  algorithm  was  also  evaluated  by  comparing  the  changes  in  total 
hemoglobin  concentration,  [HbT],  and  jus’  for  each  set  of  the  above  measurement.  Since  yeast 
did  not  induce  much  change  in  either  [HbT]  or  gs\  these  two  parameters  should  ideally  remain 
constant  when  the  oxygenation  state  was  altered.  Figure  5(c)  shows  the  measured  values  of 
[HbT]  and  ps'  (  at  750  nm  and  830  nm)  with  changing  02  Sat.  As  expected,  those  values  are 
found  to  be  constant  across  the  saturation  range  of  20-80%  with  small  SEMs  as  [HbT]  =  28.0 
±0.4  pM,  gs’  (750  nm)  =  10.46±0.05  cm"1,  and  gs’  (830  nm)  =  9.34±0.04  cm"1.  Another  set  of 
experiments  were  conducted  to  test  the  consistancy  of  gs’  when  Intralipid  volume  was  varied 
from  75  ml  to  170  ml  within  the  given  liquid  phantom  volume.  Figure  5(d)  shows  a  comparison  of 
values  obtained  by  the  ISS  oximeter  and  LRS  system  at  830  nm.  The  mean  difference  (with 
standard  deviation)  between  the  scattering  values  obtained  by  two  methods  was  calculated  to 
be  0.02±0.13  cm"1,  demonstrating  almost  perfect  match  between  the  two  methods. 

After  comprehensive  algorithm  validation  for  LRS,  we  were  able  to  quantify  optical 
properties  and  hemodynamic  parameters,  namely  ga,  gs\  [HbO],  and  [HbT],  of  interrogated 
tissues  by  the  optical  probe.  The  LRS  system  was  ready  for  studies  using  ex  vivo  human 
prostate  specimens. 
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For  AFLM  results,  we  did  not  obtain  too  much  meaningful  results,  possibly  because  our 
tissue  phantoms  are  highly  hemoglobin-based  and  are  not  a  good  light  fluorescence  tissue 
model.  As  mentioned  earlier,  the  auto-fluorescence  contrast  comes  mainly  from  Flavins,  Lipo- 
pigments  and  Porphyrins.  Our  current  tissue  phantoms  did  not  include  much  of  those 
components.  The  system  would  be  better  tested  and  validated  when  we  are  ready  for  ex  vivo 
human  prostate  specimen  measurements. 

2.3  Report  for  Aim  2 

Aim  2:  to  integrate  the  optical  fibers,  which  collect  light  scattering  and  auto-fluorescence  from 
the  prostate  tissue,  into  a  transrectal-ultrasound,  needle-biopsy  probe.  In  the  development 
phase,  the  optical  signatures  of  prostate  cancer  can  be  collected  with  the  biopsy  tissues  and 
identified  along  every  tract  of  needle  biopsies. 

For  Aim  2,  we  have  completed  and  implemented  the  1st  version  of  an  integrated  optical-fiber 
assembly  that  is  under  testing  to  collect  both  LRS  and  AFLM.  The  detailed  design  is  given 
below: 


li 


□  Fluorescence  Lifetime 

□  CW  reflectance 


700  pm  1 

Wavelength: 400  -  850  nm 

Front  view  of  the  probe  NA:  0,22 -0.4 


3.  Key  Research  Accomplishments  and  Reportable  Outcomes 

(1)  Successfully  developed  a  multi-modal,  optical  spectroscopic  instrument,  which  allows  the 
measurements  of  (1)  LRS,  (2)  AFS,  and  (3)  AFLM.  The  individual  systems  have  been  tested 
using  laboratory  tissue  phantoms. 

(2)  Two  related  research  works  are  published: 

(a)  Anindita  Mukerjee,  Rafal  Luchowski,  Amalendu  P.  Ranjan,  Sangram  Raut,  Jamboor  K. 
Vishwanatha,  Zygmunt  Gryczynski,  and  Ignacy  Gryczynski,  “Enhanced  Fluorescence  of 
Curcumin  on  Plasmonic  Platforms,”  Current  Pharmaceutical  Biotechnology,  11,  223-228 
(2010). 

(b)  Anindita  Mukerjee,  Thomas  J.  Sorensen,  Amalendu  P.  Ranjan,  Sangram  Raut,  Ignacy 
Gryczynski,  Jamboor  K.  Vishwanatha,  and  Zygmunt  Gryczynski, “Spectroscopic  Properties 
of  Curcumin:  Orientation  of  Transition  Moments,”  J.  Phys.  Chem.  B,  114,  12679-12684 
(2010). 

4.  Conclusions  and  plan  for  next  year 

In  summary,  for  the  period  from  May  2009  to  August  2010,  we  have  implemented  two  separate 
units  for  AFLM  and  LRS  and  demonstrated  that  both  AFLM  and  LRS  are  feasible  methods  to 
collect  good  optical  signals  from  tumor-simulating  tissue  phantoms.  Flowever,  the  AFLM  results 
from  tissue  phantoms  may  not  truly  reflect  the  actual  signals  collected  from  the  human 
prostates.  Further  measurements  from  human  prostate  specimens  will  be  explored  in  the 
following  year.  On  the  other  hand,  the  LRS  results  are  much  reliable  and  reflect  the  differences 
between  healthy  and  tumor  tissues  more  realistically  than  those  in  AFLM.  The  sensitivity  and 
specificity  of  this  technique  is  also  reasonably  high;  the  technique  may  allow  deeper  tissue 
measurements,  as  compared  to  AFLM,  depending  on  [HbR]  levels  in  the  tumor. 

In  the  coming  year,  we  will  carry  on  the  development  for  Aim  2  and  Aim  3:  (1)  to  integrate  the 
optical  fibers,  which  collect  light  scattering  and  auto-fluorescence  from  the  prostate  tissue,  into  a 


12 


combined,  needle-biopsy  probe;  (2)  to  collect  optical  signals  of  control  and  cancer  tissues  ex 
vivo  along  with  the  regular  human  needle  biopsy,  followed  by  classification  algorithm 
development  to  discriminate  cancer  tissues. 
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